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Abstract 



In the Starobinsky inflationary model inflation is driven by quantum cor- 
rections to the vacuum Einstein equation. We reduce the Wheeler-DeWitt 
equation corresponding to the Starobinsky model to a Schroedinger form con- 
taining time. The Schroedinger equation is solved with a Gaussian ansatz. Us- 
/\ ' ing the prescription for the normalization constant of the wavefunction given 

» ■ in our previous work, we show that the Gaussian ansatz demands Hawking 

type initial conditions for the wavefunction of the universe. The wormholes 
induce randomness in initial states suggesting a basis for time-contained de- 
scription of the Wheeler-DeWitt equation. 

I Introduction 

Inflation is an essential ingredient in modern cosmology to solve the horizon, the flat- 
ness and the monopole problems. In order to ascertain whether inflation is natural, 
there have been several attempts to study the curvature squared cosmology along 
with the Einstein curvature term. Even before Guth's [1] proposal of inflationary 
scenario, Starobinsky [2] proposed an inflationary model taking one loop quantum 
corrections to the classical Einstein Lagrangian. It has been shown that this model 
is equivalent to the curvature squared model in which the universe tunnels into the 
de Sitter phase that becomes unstable and finally emerges into the Friedmann era. 
In order to show that the inflation is quite natural phenomena for any initial 
conditions (called no hair conjecture), several authors [3] attempted to show that, 
for some models, inflation is really an attractor. The works done by Mijic, Moris 

1 



and Suen [4], Maeda [5] and Vilenkin [6-7] refer to the early universe situation 
of inflationary cosmology. Vilenkin deals with a quantum viewpoint to study the 
Starobinsky model and considers the tunneling phenomenon in the de Sitter phase 
through solution of the Wheeler-DeWitt equation [7] in quantum cosmology. In this 
approach one requires a specific ansatz for the initial conditions, that are currently 
divided into two ways: the Vilenkin proposal [8] and the Hartle-Hawking proposal 
[9]. In the approach (without phase transition), the de Sitter solution is obtained as 
self consistent solution of the vacuum Einstein equations modified by one loop quan- 
tum corrections due to quantized conformal matter fields or through the solution 
of the Wheeler-DeWitt equation with appropriate boundary conditions leading to 
an inflationary de Sitter modes. The other approach is related to phase transitions. 
Here the inflation is driven by a false vacuum energy density and is related to the 
inflaton field in the model. The no hair conjecture is found to exist in both type of 
approaches as revealed by the works of Starobinsky and Schmidt [10]. In Maeda [5] 
as well as in Mijic et al [4], the analysis of the inflation is looked upon in a classical 
way. Using Wald's [11] strong and dominant energy condition, it has been shown 
that in i? 2 -type cosmology inflation is quite natural but acts as a transient attrac- 
tor only. In Mijic, Morris and Suen the curvature squared theory has been shown 
to be equivalent to the Einstein gravity theory with a scalar field but without any 
potential form for the scalar field. Though Maeda succeeds in obtaining a potential 
form in their work with a long plateau region, it is not clear how the universe from 
a region of Planck region arrives at the plateau region simply by losing energy to 
roll over the flat region depending on initial conditions. We believe, the answer 
should come through the quantum behaviour, which is absent in their model [4,5]. 
The quantum analysis carried out by Vilenkin [7] looks nice. However the initial 
condition of tunneling behaviour with respect to no hair conjecture has not been 
tackled with justice. However Starobinsky analyzed his R 2 type theory to show that 
there is an unstable de Sitter solution followed by the present Friedman era after 
sufficient inflation. Though Vilenkin tried to answer some aspects of stability of de 
Sitter solution in the Starobinsky model, but it is not clear even now how does the 
self consistent de Sitter solution in the Starobinsky model exits to the Friedmann 
era, though it is unstable in both past as well as in the future. We try to understand 
this aspect from the energy conservation. 

Moreover the curvature fluctuation term has been introduced in an adhoc way. 
In calculating the curvature fluctuations, Vilenkin treats the wavefunction as if it 
is normalized but it is known that the normalization of Vilenkin wavefunction is 
awkward. Whereas the instability is related to the curvature fluctuations and which 
in turn is related to the randomness in the initial states it is therefore necessary 
to study the origin of randomness. As the Starobinsky solution emerges due to 
the quantum corrections, the quantum superposition principle remains inherent in 
this description. Hence the emergence of the classical de Sitter universe should 
be answered in a natural way from the study of the quantum wavefunction of the 
universe. 



In this paper we try to answer this problem at least being successful to point out 
the direction in which the investigation has to be persued to get an final answer. Our 
aim is to study (i) the normalization aspect of the wavefunction, (ii) the tunneling 
probability and (iii) the randomness and the curvature fluctuation in the Starobinsky 
description. It is shown that the energy momentum conservation forces a de Sitter 
universe to tunnel into a realistic spacetime. The tunneling from the Euclidean to the 
Lorentzian sector also demands a Wheeler-DeWitt type equation for the quantum 
wavefunction (referred as instanton) to settle the initial conditions which in turn is 
related to the two proposals of the quantum wavefunction of the universe. Here we 
find that the wormhole connected universe greatly modifies the initial conditions to 
be imposed on the Euclidean de Sitter instanton universe. We report in this paper 
the main results. 

The quantum instability of the de Sitter spacetime has been discussed at length 
by Mottola [12] and many others to understand the nature of Im L e ff with respect 
to the de Sitter spacetime. It may be pointed out that the Im L e ff is a measure of 
the instability. From these works and from the work of Nevelli [13], it is clear that 
the tunneling paths in the Euclidean sector are crucial to understand the instability 
of the de Sitter vacuum. But there are no successful attempts to include such a 
contribution in the Starobinsky type model. In section II we discuss the Starobin- 
sky model with some known results, barely needed in our follow up discussion and 
also discuss the tunneling of the universe in the de Sitter space, fixed by the energy 
momentum conservation. This treatment allows us to persue the quantum evolution 
of the universe in the context of the Wheeler-DeWitt equation and study the evolu- 
tion of the universe in R 2 cosmology. In section III we introduce a time variable to 
reduce the Wheeler-DeWitt equation to the Schroedinger form (containing time) of 
quantum mechanics. The section IV deals with the calculations of the wavefunction 
in the Starobinsky model. The section V ends with a discussion. 

II The Starobinsky model 

In the Starobinsky model inflation is driven by the quantum corrections to the 
vacuum Einstein equations 

R»u - \g^R = SnG < T^ >, (1) 

where < Tu V > is the quantum corrections and is given by 



where 



<T, U >=^ W H^ + P^H^ : (2) 



,1)„ = R;fj,;iy — 9fivR;a' U + RRfiv — ^QtivR (3) 



and 
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(3)„ = RZRva - ^RRfiu - 7:9iivR UT RuT + ~;9p,vR ■ (4) 

We use the following conventions 

Rp-ys = dsTpy — djT a ps — r a jU7 r M / 35 + r a ju5 r /l / 3 7 . (5) 

and the Ricci identity as 

A ; ^- u — A . U .^ = AK /3^ u , (7) 

with R a /3js antisymmetric in last two indices. The metric is assumed to be of the 
Robert son- Walker form, 

ds 2 =C(r})[dr] 2 -dT,l(r, $,(/))], (8) 

where dT<k is the metric on a 3-sphere, 3-plane and 3-hyperboloid for closed (k = 
+ 1) , flat (k = 0) and open (k = —1) metric respectively. As we would need some 
properties of {l) H v and (3)^ v in our description, we mention here some salient 
features. The tensor (1) H is identically conserved. Using the Bianchi identity 

(R^ - \g^R) = o, (9) 

we get from (3) 

(X)h]i;v = R ' V ;»;v ~ 9 I l [R;a;v' a + R ;v R^. (10) 

Defining K ,v = A v and using Ricci identity Eqn.(7), we get 

(!W=0. (11) 

The tensor (3) H , on the other hand is conserved only in the conformally flat 
spacetime and can be written as 

( 3 W^ = ~^ R2g ^ v ~ Rpa Rp^u- (12) 

The trace of the above two tensors is given by 

(l)^ = -3V (13) 

and 

(3)^ = \r 2 - R aP R a ". (14) 

The sum of the Eqns.(13) and (14) gives the so called trace anomaly i.e., the regular- 
ization carried out in n-dimension keeps a footprint when put back in 4-dimension 



and thus breaks the conformal invariance of the original action. In the RW space- 
time since we have a single variable C(rj), it is sufficient to consider the time-time 
components of Eq.(l); all other components will be linearly related. We write (1) H 
and (3) H explicitly. 
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(1) ==-C- 1 [-9DD + ^-D 2 + ^D 4 + 9kD 2 -l8k 2 }, (15) 
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(3) Hw =C-^D* + 3 -kD 2 + 3k% (16) 

i2oo-^ooi2=-3(^ + fc). (17) 

Here D = ^ and the dot refers to the derivative with respect to rj. Introducing 

dt 
drj 



a (ty 

we get the '00' component of Eqn.(l) as 



a' 2 + k 1 l a' 2 + k y 1 [ 2a! a!" a^_ 2a" a' 2 a' a' 2 k*_ 



a* Hq [ gt M l cr a" a° a a* a* 

(19) 
where all constants /3, oc and G are incorporated in H and M . The Eqn.(19) has 
solutions 

a(t) = Ho' 1 cosh H t, k = +1, (20) 

o(t) = a exp (if *)> fc = °» ( 21 ) 

o(t) = Ho -1 siiih f/ot, k = -1, (22) 

Eqns.(l) and (19) refer to the Starobinsky description. In quantum language it is 
said that the universe tunnels quantum mechanically from "nothing" to the de Sitter 
inflationary mode, (Eqns.(20,21,22)). Here "nothing" means no classical spacetime. 
At the moment of nucleation, t — the universe has a size a(0) = H~ l , zero velocity 
i.e., a'(0) = 0. This is the beginning of time. In quantum theory, tunneling is mainly 
due to vacuum fluctuations. The quantum cosmology deals with the tunneling 
of universe from the Euclidean region a < H~ x usually termed as the classically 
forbidden region to a classically allowed region a > H^ 1 . There are several modes of 
inflationary models, among which the de Sitter mode of evolution plays a key role. 
If the de Sitter phase plays any previledged role with signature change occurring at 
the moment of nucleation, one has to know the behaviour of quantum wavefunction 
of the universe in the Lorentzian regime from a path integral approach or from the 
tunneling approach. 



To decide some of this aspects and also the stability of the Starobinsky phase 
we note the following exercise. Let us write (1) term in two parts 

Ww v = [-^;MI" ~~ 9fJ.uR;a' a - RR^u H J ] + [2R(Rfj,v — )] 

= (!)#„, + (1)5^ (23) 

and neglect the second term to construct < T^ u >. The < T^ v > so constructed is 
not conserved as is evident from Eqn.(lO) and (11). Neglecting (1)5 v term we get 

2 Wtfoo = U -9DD + ^- - r±L + 27&L> 2 C + ISA: 2 ). (24) 

G o 2 

Using Eqn.(24) and (16) in Eqn.(l) we get for the '00' component of the Einstein 
equation as 



,/2 , j_ 1 r _/2 - -■ 2 



eT + k _ 1 



a' 2 + fc 



a 2 



M, 



2 



2a! a!" a^ 2aV 2 5a /4 a^ _ k? 

a A a z a 6 a 4 a 4 a 4 



(25) 



with a' = ^,dr) = -%:. All constants /3,a and G are incorporated in Hq and 
M . The Eqn.(25) has a characteristic feature that, it has also the same set of 
solutions (20-22) as well as the Euclidean de Sitter solution for k = +1, 0, —1 but 

< T v > 7^ 0, as is evident from the expression of (l) s v . To show it we convert 

Eqn.(25) setting H = — as 

TT 2 

H 2 {H 2 - H 2 ) = ^±.{2HH" + 10H 2 H' + H' 2 ). (26) 

M 

This shows that H = H is a solution of Eqn. (26). This result seems quite 
interesting. In realistic spacetime we know T" u = 0. So emergence in the de Sitter 
mode is possible only if (l)s = 0. From the expression given in Eqn. (23), we 
find that for the de Sitter type of solution (l)s exactly vanishes. This explains 
the emergence of de Sitter phase in the Starobinsky description. As soon as the 
universe evolves into the Starobinsky type de Sitter phase we have t — i.e., the 
universe starts at a = H~ l . As we observe the de Sitter mode of solution even with 

< T" > 7^ 0, the existence of such a solution is permitted provided the violation 
obeys 

(AT 00 )(Ar) » 1, (27) 

where At << the Planck Time. The Eqn. (27) then implies AToo >> 10 19 GeV 
i.e., we have a hot big bang scenario with temperature ~ 10 l7 GeV. 

It is quite probable that at such a short scale (10~ 33 cm) higher order corrections 
like R 3 , .R 4 etc. begin to compensate the violation. However, whatever be the 
magnitude of compensation, the violation is there and the only way the universe 
tunnels to a realistic spacetime is through the de Sitter mode of evolution making 
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(1)5 v — and there is a signature change at the moment of nucleation. The universe 
tunnels from the Euclidean de Sitter to the Lorentzian de Sitter solution. Using the 
expression for (1)5 , we get 



I_LIS 



k 2 k 2 

l AT »l K r* (28) 

taking a'(0) = around the nucleation point. Now o(0) = if -1 , so 

I AT 00 I « A; 2 # 2 . (29) 

At this point one may be tempted to argue that the tunneling to the Starobinsky 
phase occurs with higher probability when H gets smaller. But we then have to 
abandon the inflationary scenario and its fruits related to the horizon and the flatness 
problem. However there is a dynamical zeroing of |AToo| leading to de Sitter mode 
causing (1)5 v — 0. In order to incorporate this dynamical aspect in the framework 
of the quantum cosmology we lean towards the calculation of tunneling probability 
to the Starobinsky phase. The treatment will help us answer about the nature 
of the boundary conditions that would satisfy the solution of the Wheeler-DeWitt 
equation. 

It should be pointed out that while subtracting a part from the energy momentum 
tensor, as if we spoiled the energy momentum conservation still getting the de Sitter 
solution with the truncated T 00 . What we observe that only for the de Sitter solution, 
(l) s term is automatically zero and the energy conservation is restored. For any 
other solution energy momentum conservation would be violated. This is the reason 
of spontaneous nucleation of the Starobinsky phase into the de Sitter phase, hitherto 
not pointed out in any previous work. One should not confuse that we get de Sitter 
solution even violating the energy conservation. This would, if true, then leads 
to horizon size fluctuations at the onset of the FRW hot big bang destroying the 
benefits of inflation. What we like to point out that the uncertainty principle plays 
a decisive role in setting the quantum character of the universe as well as the initial 
conditions. 

Ill Time in Quantum Gravity and The Schroedinger 
Equation 

The discussion in section II reveals that the quantum corrections to the vacuum 
Einstein's equation drive the inflation in the Starobinsky model and there is a tran- 
sition from the Euclidean to the Lorentzian region. Moreover, the transition to the 
classical spacetime occurs for T v = 0, with a de Sitter inflation. There occurs thus 
a quantum to classical transition in which energy conservation plays a deciding role 
as if it acts as a boundary condition for the transition. 



In the classical description, the time and the Hamiltonian are related through 
the Hamilton's equation whereas the Schroedinger equation plays the same role in 
the quantum description. But in quantum gravity, because of the constraint relation 
H = 0, quantization of the gravitational field faces a conceptual problem in that 
no time parameter appears at a fundamental level. This is the problem of time in 
quantum gravity and recently the matter is under serious investigation by various 
authors [14-20]. 

The energy conservation, the Schroedinger equation (the Wheeler-DeWitt equa- 
tion in quantum gravity) and the interpretational framework require the introduction 
of an external time parameter to understand in a coherent way the quantum to clas- 
sical transition and the tunneling problem in the quantum gravitational description. 
We will be mainly concerned with these two aspects in the Starobinsky description. 
The emergence of universe in classical region requires the solution of the Wheeler- 
DeWitt equation and the effectiveness i.e., the interpretational framework will be 
discussed in studying the decoherence mechanism in quantum gravity. 

It has been shown by Santamato [21] how to arrive at the Schroedinger equation 
starting from classical description using Madelung-Bohm [22,23] and Feynes-Nelson 
[24,25] approaches bearing epistemological content of traditional quantum mechan- 
ics. He starts with a function S(q, t) satisfying the Hamilton- Jacobi equation of the 
classical mechanics 

QQ 

— + H(q,VS,t) = (30) 

where H(q,p = WS,t) is the classical Hamiltonian such that 

t dH(q,VS,t) 
v{q,t) = ■= (31) 

OPi 

and following Madelung-Bohm and Feynes-Nelson approaches defines a probability 
density p(q, t) satisfying the continuity equation 

d t p + d i (pv i ) = 0. (32) 

It has been then shown that the function 



ip(q,t) = y/p(q, t) exp 



h s M 



(33) 



satisfies the Schroedinger equation for a given Hamiltonian combining (30) and (32) 
with gradient of scalar curvature of space (the Weyl space) introducing randomness 
on initial positions. 

As there is no time variable in quantum gravity we demand S(q, t) = S(q) and it 
satisfies the Hamilton Jacobi equation but source free. Defining a time parameter, we 



call it WKB time, it is obvious that Jp(q, t) oc ip(q, t) then satisfies the Schroedinger 
equation when S(q,t) is time independent. Though not identical, the same view 
is also expressed by Kiefer [18,19], in defining a time operator -^ starting from a 



minisuperspace Wheeler-DeWitt equation. We discuss briefly the outline of our 
approach with minimally coupled scalar field in gravitational background. The 
details is placed elsewhere [26]. In the next section the approach will be elucidated 
for the Starobinsky description. 
We start with an action 



j4 / — 
a xy—g 



■R 



8nGJj2 



16ttG 2tt 2V 2 
d 3 xVhK 



' '\&^ + V(d)) 



in the FRW model 

ds 2 = -dt 2 + a 2 {t) 

The Hamiltonian constraint is 

1 



dr 2 



1 — kr- 



+ r 2 (d0 2 + shr 



■u,a p - + -h p * , -"' u + ' fv M = - 



The dynamical equations are 



a 

a 



and 



£ + ^ + s<^-"<*» 



h ■ dV 

+ 3-0+ — = 0. 
a ad) 



In (34) to (37), M = §, 



<->'' — Planck 



Pa 



, k — 0, ±1 and 
-Mad, Pa, = a 3 d>. 



(34) 



(35) 



(36) 



(37) 



(3? 



(39) 



The prime denotes derivative with respect to time and k is the trace of extrinsic 
curvature. Identifying P a = ^-^P^ = ^, the Hamilton- Jacobi equation is 



l^r + ^(^) 2 - \Mka' + a l V(o) = 0. 



1 0S 2 x vu 2 x „,„_;? j 

2M l da ' 2a 2 y 86 ' 2 



We define a time operator (a directional derivative) 

d_ _ ^ ( dH__d_ _ <9g <9 



(40) 



(41) 



Going to the FRW minisuperspace, we get the Wheeler-DeWitt equation Hip = 
i.e., 

h 2 , d 2 p d x ft 2 <9 2 MA;a 2 



2MW 



+ 



<9a ; 2 a 2 <9</> 2 



+ aV(0) 



^ = 0. 



(42) 



Identifying P a = —%h^ and P$ = —%h-jL and substituting 



i/)(a, 0) = exp 



oc 
n=0 



Y.AnM 1 - 



$(a,0) (43) 



in (42) we find in M 2 — order 4np = 0, and in M 1 — order 

yO 2 + ^ = ° < 44 > 

Identifying S = MA , equation (44) gives the source free Hamilton- Jacobi equation 
see (40): 

4-P a 2 + — ka A = 0. (45) 

2M 2 v ; 

Using this condition in (41) we get 

d_ = _j_ds IL l_ () 

dt Ma da da' v ; 

since S = S (a) only. Thus writing 

^(a, 0) = e lAo{a)M $(a,0) (47) 

and substituting (47) in (42) and using (44) and (46) we get 

(9$ 
ih-^-H^aA) $ (48) 

neglecting % — order term ^f . As mentioned in the beginning of this section, we 
obtain the reduction (48) using (41), instead of the continuity equation. In the next 
section we take up the Starobinsky model to evaluate \ip\ using a Gaussian ansatz 
for $(a,0). 

One might point out various drawbacks in obtaining (48) using (41), (43) and 
(46), which are now been persued by many authors in the framework: canonical 
quantization of gravity [27,28]. Recently we have been able to obtain (48) using 
the prescription of 'time before quantization' [26] . One may consider (48) from a 
different angle. Suppose we have (48) valid in the semiclassical region, we construct 
H^ for the Starobinsky description, adopt a boundary condition for (of course 
valid in the large scale factor region), investigate whether the boundary condition 
is anyhow related to the boundary condition at the small scale factor region and 
reproduces the wavefunction of the universe suiting a given boundary condition 
proposal. 
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IV Wavefunctions in the Starobinsky model. 

The main problem in effecting canonical quantization in the Starobinsky description 
in minisuperspace formalism is the absence of an action in closed form. However it 
has been shown that [7] for M 2 << H 2 an action in closed form is obtained as 

S = 2tt 2 f L(R)a 3 dt (49) 

with 

1 / „ A li , it . , 

W-utoG^wt + K^K* (50) 

where R = 12H 2 . Evaluating the action it is found that it has maximum at R — R 
so that the semiclassical tunneling probability is given by 

An 
PocexpH^I) = exp(-£^). (51) 

The curvature fluctuation is obtained by expanding S(R) around R = R 

S(R) = S + ^S"(R )(5R) 2 . (52) 

Using (49) and (50), S(R) is given by 

S{R) = ir ( R + 6Ke + M ]n R, ) (53) 



so that the curvature fluctuation 

, 5R . n GH 



Ro 7T 



(54) 



It is interesting to note that this curvature fluctuation is proportional to |AT 00 | 
as is found in (28). However the result (54) is classical. An informative description 
for the nucleation of the universe can be obtained by solving the Wheeler-DeWitt 
equation for the wavefunction of the universe. Vilenkin [] attempted this problem 
solving (42) for the Starobinsky model with the tunneling boundary condition to 
obtain a curvature fluctuation of the order of (54). Our aim is to study these aspects 
with the formalism just mentioned above to explore: 
(i). the normalization aspect of the wavefunction 
(ii). the tunneling probability 
(iii). the randomness and the curvature fluctuation 
(iv). the randomness and the decoherence and 

(v). the backreaction and the wormholes contribution to the solution of the Wheeler- 
DeWitt equation. In this paper we take up the first three points. 
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Using the method outlined in [7], the Wheeler-DeWitt equation reads 

d 2 1 d 2 



dq 2 q 2 dx 2 



-V(q,x) 



^(q,x) = 



(55) 



where q = H a 



v \k 



with V = L'(R) and L ' = L'(R ), and 

X = 2 l ^R~^ 
V(q,x) = X~ 2 q 2 (l-q 2 +fl 2 ( x )q 2 ), 

M 2 
H*{x) = -^(2x + e- 2x -1), 



and A 



2H 

GM 2 
6tt 



(56) 

(57) 
(58) 

(59) 



M n4 



In deriving (55) it is assumed that M « H and terms ~ {j^Y are neglected. 
Assuming \x\ < 1 and substituting Q = -4=, we find using (55) and (57) 





" d 2 

OQ 2 


-QW-^-^H 


where 




U{x) = A(l - mV) 


and 




2 M 2 

m — 



* = 



H. 



2 ' 



(60) 



(61) 



(62) 



As discussed in the previous section, we convert (55) in the Schroedinger form 
substituting the WKB form 



^ = exp(^M 1 - n A„)^(g,x). 



(63) 



Collecting terms in different order in M we obtain for M 2 — order ^tt 9 - = 0, which 



implies that A is purely a functional of the gravitational field and M 1 — order gives 
the source free Hamilton- Jacobi equation 



dQ' 



We introduce a time operator as 



d 



Q 2 = 0. 



1 0A o d 



dt Q dQ dQ 



(64) 



(65) 



and substitute (63) in (60) to get 
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,dip(Q,x) 

%- 



dt 

r e 1 
have set M = 1. Thus we have 



1 d 2 Q\ u . 1 



2Q 3 dx 2 2 w 2Q 

d 2 t/ 



1>{Q,x). (66) 



In obtaining (66) we have neglected terms ~ -^. After the WKB reduction we 

(67) 
(68) 
(69) 



^ o ( q ) = e ± 2x. (70) 

The solution with the negative sign in the exponent was obtained by Vilenkin with 
the assumption that q << 1 and 

— (0,x)=0 (71) 

and the tunneling solution exponentially decreases in the classically forbidden region. 
From (55) the classically forbidden region lies between q = and q ~ 1 + t>/J> 2 (x) 





q = MQ)*P(Q,x 


where 


^ o = e *MQ). 


From (64) and (68) one finds 






MQ) = e ± ^. 


In terms of the q variable 





and acts as turning points. 

Now to solve (66) we make a Gaussian ansatz 

^ = N(t)e- i2 ^ x2 . (72) 

Substituting this in (66) leads to coupled equations for Q and N 

<<„,*= n +0 . A + £ (73) 



*>= "'* ■ (74) 



and 

tt 2 + Q 6 \m 2 

~Q~ 3 

In (73) and (74), t is the WKB time defined by (65) and parameterizes the classical 
trajectories in the minisuperspace which is now spanned by Q. With the ansatz 



^y 



one finds from (74) 



n = -iQ 6 - (75) 

y 



V + 3 qV ~~ Xm2y = °' ^ 
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Introducing the conformal time n by the relation dt = Qdr), we get from (76) 



V" + 2^2/ - \m 2 Q 2 y = 



(77) 



where prime denotes a derivative with respect to n. We now specify the model by 
specifying Q = — j=- such that the background undergoes an exponential expansion. 
Equation (76) now reduces to 



2 , m 2 

-y — jy = ° 



and is solved by 



y = v 



3/2±y / 9/4H 



(78) 



(79) 



Approximating J9/4 + m? = (3/2 + m 2 /3) (since m 2 = jf? « 1), and using (75) 
in conformal time coordinate we get, 



n 



-im 2 v\ 



3 ' 



(80) 



As Q is imaginary the state (72) will not be normalizable. One of the ways to obtain 
the real part in Q is to consider various mode solutions of the scalar field x as in the 
work of Kiefer [3] . However we will follow a different line based on our previous work 
[29] having interpretational significance along with settling the boundary conditions. 
Remembering that m 2 y\ is a very small quantity, we approximate (73) for large Q 

i^-lnN = tt 3 \. (81) 

dt 



With | = V\Q4j, one finds from (81) 



dQ 



N = N exp 



iQ 3 V\ 



so that 



■ip = N exp 



iQ 3 V\ im 2 V\Q 3 



6 



-x 



33) 



where N is a constant to be determined. The bracketed term in the exponent (83) 
is approximated as 



S, 



eff 



iQ 3 V\ irrV 



v/AQ 3 



6 
Q 2 \(l-fi 2 (x)) 

Q 2 \(l-fi 2 )-l 
14 



-x 



3/2 1 + /i 2 



A 

3/2 1 +/i 2 (x) 

X 



14) 



In our previous work [29] we have shown that the multiple reflections between the 
turning points will contribute to the normalization factor and are supposed to arise 
from the wormhole contributions. Using (84) and the result of [29] the normalization 
constant N , according to wormhole dominance proposal, turns out to be 



N n 



expS eff (Q x ,0) 
l-exp[2S e// (Q,0)] 



where 



S e ff(Qx,0) - S e ff(Q)\ X 



(85) 



(86) 



and Q = and Q = Q T = — 

(A(1-M 2 M))5 

(55) and (57). We thus get from (85) and (83) 



1 



- are the turning points as can be seen from 



< ' - ?w {+3^(1 + A*)) [i - Kq - </V - i) 3/2 J } 



(87) 



For simplicity we have dropped the denominator of (85) but is nonetheless impor- 
tant while discussing the normalization. Eqn.(87) continued in classically forbidden 
region gives 



i> 



ex p{js( 1 + / i2 ( a 0) 



1 



l-r + q l \i 



2..2N3/2 



} 



(l-exp{2/3A(l + /i 2 (x))}) v ' 

Eqn.(87) and (88) are the wavefunctions that one gets from the Hawking's proposal. 
We can now discuss the curvature fluctuations in the newly-born universe. From 
(87), the probability of nucleation with a certain value of x is proportional to 



oc exp 



2fi 2 {x) 
3A 



(89) 



In obtaining (89), a factor (e 
(89). For x « l,fi 2 {x] 



:(1+M 2 ) 



2 2 

m x 



■M Z 



oc exp 



l) 2 in the denominator is kept as multiplying 
)x 2 and we get 



Attx 2 
GH A 



Now it follows from (56) that x 
fluctuations. Hence we can write 

,8R~ 



|g, where SR 



< i^) >=4<x 2 >-- 

tin 



GH 

2tt 



(90) 
R — R Q is the curvature 

(91) 



Exactly this result was obtained by Vilenkin from the tunneling proposal, but our 
wavefunction corresponds to the Hawking's proposal. 

Though the result (91) is identical with the Vilenkin, there are some marked 
interesting differences. In Vilenkin [7] tjt(1 + H 2 (%)) term is added in (87) to get 
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the behaviour S — > %- for q — ► 0. In our approach when the denominator in (87) 

2 

is taken into account, S — > — |r as g — »■ 0. If we set aside the denominator in 
(87), the wavefunction looks apparently as the Vilenkin wavefunction. By analytic 
continuation, in the classically allowed region, Vilenkin gets back the term — ^(1 + 
/i 2 ) in the exponent to find 

V> ~ exp {--L(l + fi 2 (x)) [l - t(q 2 - oV - l) 3 / 2 ] }. (92) 

This analytic continuation seems questionable. Thus we see that the Gaussian 
ansatz and the normalization prescription serve as the boundary conditions to settle 
between the two current proposals in favour of the Hawking's prescription. We 
also stressed this in our previous work [29]. The present work works out explicitly 
the equivalence between the timeless Wheeler-DeWitt equation (55) and the time 
contained Wheeler-DeWitt equation (66) and hence the time parameter prescription 
as mentioned in this work suits also in the frame work of the quantum cosmology. 



V Discussion 

We find that the normalization factor arises from the repeated reflections between 
the turning points and hence corresponds to the higher order corrections related to 
the quantum fluctuations in the WKB description. In our previous work [29] we have 
shown that the factor N can be interpreted as contributions from the wormholes 
using Klebanov and Coleman's arguments [30,31]. Apart from the wormhole picture, 
it can be said that the quantum force has its origin in the curvature fluctuations 
with 'to and fro' motion within < q < 1 + |/x 2 (a;) and it necessitates to uphold the 
probabilistic interpretation. Leaving aside the interpretational hindrance, 'to and 
fro motion' is assigned to terms like e lS and e~ tS but the Gaussian ansatz required 
only the form (72) in the classically allowed region as if there is a suppression of 
interference terms. It has been argued that some boundary conditions at small scale 
[32] would lead to the quantum effects in the vicinity of the turning points. Thus 
the only way to interpret e tS and e~ tS superposition [i.e., the quantum effects] is 
to turn back towards Klebanov and Coleman's arguments of wormholes connections 
around q ~ 0. There is also an objection about the 'to and fro' motion with re- 
spect to time because there is no external time parameter with respect to which 
the universe can turn. But our calculation in the time variable description correctly 
obtains the wavefunction both in the classically allowed and in the forbidden re- 
gion. The classical turnings points serve as a clue to obtain the peculiar behaviour 
of wave packets at a late time (q >> 1) whereas the quantum turning point (we 
identify it to be q = 0) leads to 'to and fro' quantum fluctuations. The Euclidean 
and the Lorentzian time both are manifestation of the spacetime structure latent 
in the description of the Wheeler-DeWitt equation. There is also a question about 
the influence of excited matter states and whether they lead to decoherence. In 
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the present work we prove the decoherence in reverse way i.e., leading to the form 
(88) from the Gaussian ansatz (72); however we have been able to show [33] that 
the presence of effective decoherence even in the Starobinsky description, using the 
technique of ref.[16]. Perhaps the arrow of time starts functioning from the classi- 
cal turning point when the universe emerges from the Euclidean to the Lorentzian 
regime with a signature change. It should be pointed out that the denominator 
in (88) arises from the multiple reflections from the turning points and indicates 
decisively in favour of the Hawking wavef unctions. If we accept the Klebanov and 
Coleman's [30,31] arguments, it would not be unjustified to comment that worm- 
holes contribution generates the initial randomness, and the quantum behaviour 
thus generated even persists in the classical spacetime. The emergence of the classi- 
cal universe is then explained through the decoherence to suppress the interference 
of expiS* and exp—iS like terms. The code of the decoherence is seeded in the 
normalization constant rather than in the initial conditions. 

It is worthwhile to point out that, at the semiclassical level, the Starobinsky infla- 
tionary scenario has been criticized [34] on the ground that the inflationary solutions 
are not perturbatively expandable in the parameter of "quantum corrections" terms 
in the equation of motion or the Lagrangian. There are also some questions of treat- 
ing the higher order terms in the Lagrangian on an equal footing with the Einstein 
terms. Even then, studying the quantum cosmology of the Starobinsky description 
from the standpoint of quantum to classical transition would help us rethink about 
the criticism labelled against the model. Leaving aside this fact, it is instructive to 
look at the model as a toy example, to understand the current boundary condition 
proposals, decoherence mechanism and also the origin of the quantum force in the 
early universe. 

Further there are also some drawbacks in expanding the Wheeler-DeWitt wave- 
function as a power of Planck mass as in (43) and in obtaining the Schroedinger- 
Wheeler-DeWitt equation (SWD) from it. Recently [26,35] we have been able to 
obtain the SWD equation without using the WD equation and also the expansion 
(43) and it has been found that the SWD wavefunction can be defined on the stan- 
dard Hilbert space of quantum mechanics. Though the Starobinsky description itself 
has some drawbacks as discussed in [30], our attempt is to understand the quan- 
tum to classical transition in the framework of the model and to investigate which 
of the boundary conditions proposal suits to the inflationary description. We find 
that (i) the classical Starobinsky model is consistent with the quantum description, 
(ii) the wormhole dominance proposal [29] correctly connects the Wheeler-DeWitt 
wavefunction and the Schroedinger- Wheeler-DeWitt wavefunction in the respective 
regime, (iii) the curvature fluctuation is correctly reproduced from the wormhole 
dominance wavefunction, (iv) the wormholes initiate the quantum randomness in 
the initial stage (i.e., < a < H~ v ) and (v) the decoherence is effectively [33] repro- 
duced in the Starobinsky description, provided we accept the wormhole dominance 
proposal keeping intact the fruits of the inflationary scenario. 
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